Extended quintessence and the primordial helium abundance 



Xuelei Chen* 1 , Robert J. Scherrer^ 1 ' 2 , Gary Steigman-t 1 ' 2 
1 Department of Physics, The Ohio State University, Columbus, OH 43210 
2 Department of Astronomy, The Ohio State University, Columbus, OH 43210 
(February 5, 2008) 

In extended quintessence models, a scalar field which couples to the curvature scalar R provides most of the 
energy density of the universe. We point out that such models can also lead naturally to a decrease in the 
primordial abundance of helium-4, relieving the tension which currently exists between the primordial helium-4 
abundance inferred from observations and the amount predicted by Standard Big Bang Nucleosynthesis (SBBN) 
corresponding to the observed deuterium abundance. Using negative power-law potentials for the quintessence 
field, we determine the range of model parameters which can lead to an interesting reduction in the helium-4 
abundance, and we show that it overlaps with the region allowed by other constraints on extended quintessence 
models. 



I. INTRODUCTION 

A great deal of observational evidence currently points to- 
ward a cosmological model with a nonzero cosmological con- 
stant A. A combination of the supernova la measurements 
[1], measurements of the baryon fraction in galaxy clusters 
[2], and the location of the first acoustic peak of the mi- 
crowave background anisotropies suggests a model with 
n M = 0.3 - 0.4 and fl A = 0.6 - 0.7. 

A genuine cosmological constant is a disaster from the 
standpoint of particle physics, in which the most natural value 
of A is many orders of magnitude larger than the "observed" 
value (see Ref. [Qj for a review). However, even if a plau- 
sible particle physics mechanism were developed to produce 
such a small value of A, a second problem remains: why are 
J7a and 51m comparable today? Since J1a and Hm scale very 
differently with the cosmological expansion factor, this coin- 
cidence suggests that we live in a very special epoch. This 
problem has been dubbed the "cosmological constant coinci- 
dence problem" to distinguish it from the more fundamental 
cosmological constant problem [S,^]. 

A possible solution to the coincidence problem is to assume 
that the apparent cosmological constant is not, in fact, a true 
constant vacuum energy density, but instead is due to the en- 
ergy density of a scalar field <j>, a possibility which has come 
to be known as quintessence. It might then be possible to cou- 
ple the behavior of this field to the background matter den- 
sity in such way as to achieve 'naturally' the desired result, 
namely f2^ sa Qm today. Of course, in this case, the scalar 
field may have an equation of state intermediate between mat- 
ter and a cosmological constant. This possibility was explored 
in some detail by Zlatev et al. |^] (see also reference [^]) who 
argued that a certain class of solutions ("tracking solutions") 
will evolve toward the desired behavior independent of the ini- 
tial conditions. Although such models still require fine tuning 
to produce w Qm today [Q], they seem to be pointing in a 
more plausible direction toward this result. 

These investigations assumed minimally-coupled scalar 
fields, but it was soon realized that coupling the scalar field 



to the curvature scalar R opens up another range of possibili- 
ties, dubbed "extended quintessence." Uzan [g| examined the 
general case of a non-minimally coupled scalar field evolv- 
ing in an exponential or power-law potential, and Amendola 
explored a general class of couplings and potentials. Non- 
minimally coupled models suffer from the potential problem 
that the gravitational constant G varies in time ||[l0]], but they 
have nonetheless received a great deal of recent attention as 
possible models for quintessence Jll| ] - jisfa. 

In this paper, we point out an interesting consequence of 
non-minimally coupled quintessence models: under some cir- 
cumstances these models can lead to a reduction in the primor- 
dial helium abundance. This is of interest because of a "ten- 
sion" between the recent estimate of the primordial deuterium 
abundance, D/H = 3.0 ± 0.4 x 10" 5 [Q, which leads to an 
estimate for the baryon/photon ratio of -q = 5.6±0.5x 10 °, 
and the corresponding SBBN predicted primordial helium-4 
abundance of Yp = 0.248 ± 0.001. In contrast, the actual pri- 
mordial helium abundance is likely lower. For example, Olive 
and Steigman found [fX3|] 



Y P = 0.234 ± 0.003(stat.), (1) 
while Izotov and Thuan obtained a higher value Jlr^ ] 

Y P = 0.244 ± 0.002(stat.). (2) 

While this apparent discrepancy is insufficient to discard 
SBBN, it certainly represents a "tension" in the model. Fur- 
thermore, in addition to the three standard model neutrinos, a 
sterile neutrino may be needed to explain the results from neu- 
trino oscillation experiments [17|. If either or both of these 
were confirmed, or if there are any other light particles in the 
Universe, the breach between theory and observation on 4 He 
could become even wider. 

In light of these SBBN results, any "natural" mechanism 
which lowers the primordial helium-4 abundance must be re- 
garded as interesting. In this paper we show that extended 
quintessence models provide just such a mechanism. In Ref. 
Jll]], Perrotta et al. investigated the BBN constraint on the 
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extended quintessence model. They showed that the model is 
not ruled out by increasing helium, but they did not consider 
the more interesting possibility of decreasing helium. 

In the next section, we review the evolution of the scalar 
field in these models and calculate the reduction in the abun- 
dance of helium-4 as a function of the model parameters. 
We determine whether an interesting reduction is consistent 
with other constraints on extended quintessence. Our results 
are summarized in Sec. 3. We shall use units in which 
Ti = c = 8ttGn = 1 throughout, unless otherwise stated. 



II. THE NON-MINIMALLY COUPLED QUINTESSENCE 
MODEL 

The action for an extended quintessence model is given by 

m 



1 



f(R, Q) - -w(Q)Q i ' 1 Q iM - V{Q) + L 



fluid 



(3) 



where R is the curvature scalar of the space-time and Q is a 
scalar field. We adopt the Non-Minimally Coupled (NMC) 
model of Refs. [00, in which lo{Q) = 1, f(Q,R) = 
\F(Q)R, and F{Q) takes the form, 



F(Q) = i-aQ 2 -Ql), 



(4) 



where Qo is the value of Q today. (The zero subscript will 
refer throughout to quantities evaluated at the present). Eq. 
(0) ensures that F(Q) = 1 today. 

The coupling constant £ in the NMC quintessence model 
is constrained by solar system limits on ujjbd, the Jordan- 
Brans-Dicke parameter |ffT 



F Q 



{dF/dQ)l 



> 500, 



(5) 



as well as the experimental limit on variation of the gravita- 



tional constant pi 1 ] 



; xpe 
Q 



\G/G\ = \F/F\ < 10- u yr _1 - 



(6) 



Using the form for F(Q) from Eq. 
into: 



h, these limits translate 



£Qo < 0.022, (7) 
from the solar system constraint, and 

2£QoQo < l0- n yr-\ (8) 

from the limits on the time variation of G. 

In this model, the expansion rate of the Universe is 



h2 = Jf ( Pf + \® 2 + V{Q) + 6HCQ ® 



where pf = p m + p r is the density contribution from mat- 
ter and radiation (including neutrinos). The evolution of Q is 
governed by 



Q + 'IHQ + |^ + £RQ = 0. 
oQ 



(10) 



The scalar curvature is given by 

R = 6(H + 2H 2 ) 

= i(pf- ZPf -Q 2 + 4V + 6£(Q 2 + QQ + 3HQQ) 



As a specific case, we will consider the inverse power law 
potential for Q: 



V(Q) - V a Q- a . 



(12) 



However, our final results can be generalized to other forms 
for the potential. It is known that with this type of potential, 
' the quintessence energy redshifts more slowly than the radia- 
tion in the radiation-dominated era, and more slowly than the 
matter in the matter-dominated era, ultimately dominating at 
late times [p|~|7| Jl3| , p^| |, thus providing an explanation for the 
observed accelerating expansion of the Universe The vac- 
uum energy observed today is determined by Vq but is not 
(or is only weakly) dependent on the initial value of Q, thus, 
partly addressing the "fine tuning problem" [pUq]. 
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FIG. 1. The evolution of Q for a potential V = VoQ 4 , where Vo 
is chosen to satisfy the boundary conditions given in equations (|l4|), 
and £ = 0.007. 



(9) 

A viable NMC cosmological model must satisfy the follow- 
ing boundary conditions at a = 1: 
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FIG. 2. As in Fig. [jj for the evolution of F. 
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FIG. 3. As in Fig. Ml for the evolution of the fluid and Q-field 
energy densities. 



• The effective gravitational constant today must be equal 
to the measured value, which translates to the require- 
ment 

^0=1 = 1. (13) 

This condition is satisfied automatically by our defini- 
tion in equation (^). 

• After evolving the model from an early epoch to the 
present one, we must recover the present value of cos- 
mological parameters, which we take to be: 

Q, Ma = 0.3, H = 65 km/sec/Mpc (14) 

where, as noted above, we take fio = 1. As an example, 
we show the evolution of Q, F, and the densities of interest 
which satisfy our boundary conditions in Figs. 1-3, for the 
case V = VoQ~ A , £ = 0.007. 

The general behavior of the NMC quintessence model with 
this potential was discussed in Ref. [jo]; here we review it 
briefly. At sufficiently early times, the £RQ term in equa- 
tion dominates the dV/dQ term, and the field Q settles 
down to a slow roll regime, with an effective potential (,RQ 
balanced by a "frictional force" —3HQ. (This dominance of 
the £i?Q term has been dubbed the "R-boost" [[13]]). In this 
regime, Q is nearly constant (as seen in Fig. |l]), and the poten- 
tial term V is sub-dominant; V affects neither the expansion 
nor the Q-evolution significantly. The slow roll holds until V 
becomes significant, after which Q begins to roll fast, and the 
quintessence kinetic and potential energy come to dominate. 



In this regime, the £i?Q term can be neglected, and the field 
behaves as a minimally coupled "tracker" field with a negative 
power law potential [^[-|7| Jl^ , p^[ ] . The results we show in Figs. 
1-3 are in agreement with those of Ref. [|l3|]. 

Depending on the choice of the initial conditions, the his- 
tory of the Q field may have a few more twists than described 
above; we refer the reader to Ref. [|l3|] for more details. For 
our purpose, it is sufficient to note that the early slow roll 
and the late tracker behavior are common features for most 
of these models, and we will confine our investigation to the 
range of parameters for which this behavior holds. 

Our model appears to have four free parameters: Vb, a, £, 
and Qi (the initial value of Q), which together determine Qo. 
Note, however, that Qo is almost completely independent of 
Qi, as a consequence of the tracker behavior of the model 
This is displayed in Fig. |] 

Hence, Qo is effectively only a function of a, Vb, and £. 
However, once we fix a, and $, the value of Vb is fixed by 
boundary conditions of Eq. ([14]). So effectively, a and £, 
along with the boundary conditions, fix Qo- This dependence 
is shown in Fig. || It is obvious from this figure that Qo is 
also nearly independent of £. This follows from the fact that 
the term containing £ in Eq. (|o|) is subdominant once tracking 
behavior starts. Hence, it is a good approximation to take Qo 
to be a function of a alone for boundary conditions fixed as in 
Eq. 

Since the solar system limit (Eq. [7]) is a constraint on £ and 
Qo, it translates into a constraint on £ for a given value of a. 
This limit is shown in Fig. |] For the models we consider 
here, the limit from equation ^ is subdominant and can be 
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ignored. 

Now consider the effect on the expansion rate. Let H be the 
Hubble parameter in our NMC quintessence model, and let H 
be the Hubble parameter in a fiducial (no quintessence) model, 
which we take to be a standard ACDM model with fl\ = 0.7, 
fio = 1. (Our specific choice of fiducial model is relevant only 
at late times and does not affect our BBN calculations, but 
we choose this particular model for definiteness). In order to 
reduce the helium abundance, we would like H < H, which 
implies F > 1, during BBN. However, when V dominates, Q 
rolls "down hill" in the positive direction (assume Q > 0), so 
Qo > Qbbn- Therefore, we find F > 1 if and only if £ > 0. 

The change in expansion rate can be parametrized by a 
speed-up factor £: 



(15) 



The evolution of £ — 1 in the model of Figs. 1-3 is shown 
in Fig. [7|. The region at the right side of the graph for which 
( — 1 > arises after the quintessence field enters the tracker 
regime; in this regime the scaling of pq with a is different 
than the scaling of pa in our fiducial model. This behavior is, 
however, irrelevant for BBN; the value of ( — 1 during BBN 
depends only on the radiation content of our fiducial model. 

During BBN, p f dominates the density in equation (^), so 
the speed-up factor is determined primarily by the effective 
gravitational constant F. Since the Q field is almost frozen 
during the radiation dominated era, we have Q m Qi, so 
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FIG. 6. The strongest coupling for a given a allowed by solar 
system experiments. 
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FIG. 7. As in Fig. 1, for the evolution of the speed-up factor £. 
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FIG. 8. The value of £ — 1 during nucleosynthesis as 
a function of a; the five curves are, from top to bottom, 
£ = 0.002, 0.004, 0.008, 0.012, 0.016, 0.02. 



for (-1<1. We are primarily interested in the case Qi << 
Qo, because it is in this regime one expect a large negative 
( — 1. In this case, we get 



(17) 



As we have seen, Qo is effectively only a function of a, with 
a slight residual dependence on £ for the boundary conditions 
fixed in Eq. <Jl~4|) . Hence, the speed-up factor will be a function 
of £ and Qo alone; we expect ( — 1 oc £ from Eq. (|l7|), and 
Fig. |] shows that Qo, and hence the magnitude of £ — 1 should 
increase with increasing a. The dependence of £ — 1 as a 
function of £ and a is shown in Fig. g. 

Now consider the effect that the change in H will have on 
the primordial helium abundance. For all values of a consid- 
ered here (a < 15), F is essentially constant during BBN, 
so that we can take £ — 1 to be constant. It is then straight- 
forward to calculate the change in the predicted primordial 
helium abundance (compared to SBBN). For small changes in 
the expansion rate at BBN virtually all the neutrons available 
when BBN begins are incorporated in helium-4, so the helium 
abundance is directly related to the neutron abundance. The 
faster the expansion the more neutrons are available and the 
more helium synthesized. A slower expansion has the oppo- 
site effect. For small deviations from the SBBN expansion 
rate 

Ay w 0.08(C 2 - 1) « 0.16(C - 1) -0.08fQ§. (18) 

In general, AY is a function of both a and £. However, 
we would like to find the largest possible value of AY that 



is consistent with the solar system bounds of Eq. (j7|). Since 
AY increases with £, we use the largest possible £ for a given 
a (shown in Fig. ||) to determine AY for a given a. This is 
shown in Fig. |^. 

It is obvious from this figure that — AY is an increasing 
function of a. For instance, with a = 10, we find AY w 
—0.01, corresponding to £ = 0.004. Often the speed up in 
the expansion rate is parameterized in terms of an equivalent 
number of "extra" neutrinos, AN„. 



AN„ = ^(e-l)*-ftQl 



(19) 



This reduction in Yp from its SBBN value corresponds to 
a reduction of N v from its standard model value of 3 by 
AJV„ ps 0.74. There is another, subdominant effect of a 
slower expansion in that there will be more time available 
to destroy deuterium as well as to synthesize 7 Be which will 
later electron capture and add to the abundance of 7 Li. As a 
result, the same deuterium and lithium abundances will cor- 
respond to a slightly smaller baryon-to-photon ratio which, in 
turn, will yield a slightly smaller predicted helium abundance. 



III. DISCUSSION AND CONCLUSION 

We have investigated the effect of NMC quintessence 
model on BBN. In models with £ > 0, the expansion rate 
is smaller during BBN, and less helium is produced. The 
amount of this reduction depends primarily on the coupling 
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constant and the slope of the potential, and also very weakly 
on the initial value of the field. However, the coupling con- 
stant £ is limited by solar system experiments. We find that 
for negative power-law potentials of the form V = VoQ~ a , 
the magnitude of the reduction increases with a. This is not 
a desirable state of affairs from the point of view of reduc- 
ing primordial helium. The reason is that negative power-law 
potentials give an equation of state w — p/p in which w ap- 
proaches as a increases. In Figure 10 is shown the rela- 
tion between w and a. For instance, a — 6 corresponds to 
w = —0.43, and a variety of observations argue against a 
value of w much larger than this (see, e.g., Ref. jl9|]). If we 
take a < 6 as a reasonable upper bound on such models, then 
we find AY < 0.0061. This is an interesting reduction in pri- 
mordial helium (from Y P = 0.248 to Y P = 0.242), albeit it 
requires us to push all of the parameters in the model to the 
extreme acceptable limits. For example, while a reduction as 
small as AY w 0.002 — 0.004 would be sufficient to reconcile 
the O'Meara et al. deuterium abundance [ 14 1 with the Izotov 



and Thuan helium value fll6|], a much larger reduction would 
be required if the Olive and Steigman [ 15 1 helium abundance 
were adopted. 

If there is indeed a breach between the observed helium 
abundance and the predictions of SBBN it could be healed by 
this NMC model. Alternatively, a reduction in Y P of 0.0061 
corresponds to a shift in the equivalent number of neutrinos 
by AN V = -0.47. 

We have shown that certain classes of NMC quintessence 
models lead to a "natural" reduction in the primordial helium 
abundance. Although we have concentrated on the negative 



power law potentials, our results generalize easily to other po- 
tentials. In particular, one of the major problems with expo- 
nential potentials is that they lead to an overproduction of he- 
lium [|T|. This problem could be ameliorated in the way we 
have suggested. Of course, the tracker solution for such mod- 
els genetically predicts w — at present [|T]] and thus these 
models have problems matching other observations. We also 
expect that these results would be changed if we used a differ- 
ent coupling between Q and R than that in Eq. ([I]). Thus, it 
might be possible to obtain even larger values for AY in such 
models. 
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